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We show that it is possible to define a Lorentz-covariant reduced spin density matrix for massive 
particles. Such a matrix allows one to calculate the mean values of observables connected with spin 
measurements (average polarizations). Moreover, it contains not only information about polarization 
of the particle but also information about its average kinematical state. We also use our formalism 
to calculate the correlation function in the Einstein-Podolsky-Rosen-Bohm type experiment with 
massive relativistic particles. 
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I. INTRODUCTION 

Relativistic aspects of quantum mechanics have re- 
cently attracted much attention, especially in the con- 
text of the theory of quantum information. One of the 
important questions in this context is how to define the 
reduced spin density matrix. Such a matrix should en- 
able one to make statistical predictions for the outcomes 
of ideal spin measurements which are not influenced by 
the particle momentum. We consider this problem in 
detail in the case of massive particles. The reduced spin 
density matrix is usually defined by the following formula 

Q: 

T a x = J dp(k) (k,m,s,a\p\k,m,sX) , (1) 

where p denotes the complete density matrix of a 
single particle with mass m, dp(k) = |m is the 
Lorentz-invariant measure on the mass shell and four- 
momentum eigenvectors \k, m, s, A) (i.e., P M \k, m, s, A) 
■ fc, to, s, A)) span the space of the irreducible represen- 
tation of the Poincare group. They are normalized as 
follows 

(p, to, s, a\k, to, s, A) = 2fc°(S 3 (k — Tp>)5 a \. (2) 

The action of the Lorentz transformation A on the vector 
\k, to, s, A) is of the form 

[7(A) \k, to, s, A) = V s aX (R(A, k)) \Ak, to, s, a) , (3) 

where T> s is the matrix spin s representation of the SO(3) 
group, R(A, k) — Ljj^ALk is the Wigner rotation and 
Lfc designates the standard Lorentz boost defined by the 
relations L^k = k, Lz — I, k — (to, 0). 

The key question is whether the reduced density ma- 
trix is covariant. In 11 it was stressed that the matrix 



(|T|l is not covariant under Lorentz boosts. It means that 
when we calculate the complete density matrix as seen 
by the boosted observer 

p' = U(A)pUHA) (4) 

and then the reduced spin density matrix t' oX (using Eq. 
(|T|) with p replaced with p') we find that we cannot ex- 
press t' only in terms of r and A. The reason is quite 
obvious — the Wigner rotation in the transformation law 
© is momentum dependent, except of the case A 6 0(3). 
From the group theoretical point of view it is related to 
the fact that the Lorentz group and the rotation group 
are not homomorphic. Notice that in the nonrelativistic 
quantum mechanics it is possible to define the Galilean- 
covariant reduced density matrix by the formula analo- 
gous to Eq. 2] because such a homomorphism exists. 



II. COVARIANT REDUCED DENSITY MATRIX 

As was pointed out in Q matrix is not always rel- 
evant to the discussion of relativistic aspects of polar- 
ization experiments (see, however, Q). For this reason 
we propose here another definition of the reduced density 
matrix. This definition relies on the analogy with the po- 
larization tensors formalism used in quantum field theory. 
As a result we obtain the finite-dimensional matrix which 
contains not only the information about the polarization 
of the particle but also the information about average 
values of the kinematical degrees of freedom. Moreover, 
such a matrix transforms covariantly under the Lorentz 
group action. 

To begin with we introduce vectors |a, k) such that 

\a, k) = v aa {k) \k,m, s,a) (5) 

which are assumed to transform under Lorentz transfor- 
mation A due to the following, manifestly covariant, rule 
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where D(A) is a given finite-dimensional Lorentz group 
representation. Consistency of the rules EJ) leads 
to the Weinberg-like condition [jj, ||| which has to be 
fulfilled: 

D(A)v(k)V sT (R(A, k)) = v(Ak), (7) 

where v(k) denotes matrix [iw(/c)]. Thus to calculate 
v{k) it is enough to determine v(k) and use the formula 
v(k) = D(Lk)k which is a consequence of Eq. 0. 

Assuming that the condition can be solved we can 
define the following (unnormalized) covariant reduced 
density matrix: 

a p = / dn(k)(0,k\fi\a,k). (8) 

We can easily check that this matrix is manifestly covari- 
ant under the transformation 0, namely we have 

9' = D(A)6»D t (A), (9) 

where 9 — \9 a p\. 

One can also easily verify that the matrix is Hermi- 
tian and positive scmidefinitc (similarly as 0). Trans- 
formation preserves Hermicity and positive semidefi- 
niteness of 9 but changes its trace. 

It is clear that we can define also normalized density 
matrix 

'~4i- < 10 > 

Such a matrix transforms according to the rule 

9 = D( f Dt(A . (ii) 

Tr(6»D t (A)D(A)) 

One can check immediately that Eq. I|ll[) gives a nonlin- 
ear realization of the Lorentz group connected with the 
quotient space SO(l, 3)o/SO(3). Therefore this realiza- 
tion is linear on the rotation group. However, to extract 
information about polarization of the particle it does not 
matter which matrix we use, 9 or 9. Moreover, when we 
consider representations of the full Lorentz group (i.e., in- 
cluding inversions) the most convenient choice is to con- 
sider the matrix 

n = 9T, (12) 

where L fulfills the condition 

D^T = TD -1 , (13) 

which means that in this representation P represents 
space inversions. Thus the matrix f2 transforms under 
the Lorentz group action in the following way: 

Q' = D(A)ftD~ 1 (A). (14) 

We see that transformation i|14|) does not change the 
trace of Q. Of course, having Q we can easily determine 
9 and normalized density matrix 9. 



Hereafter we restrict ourselves to the case of a spin- 1/2 
particle; generalization to the higher spin is immediate. 
In this case the Weinberg condition can be easily 
solved. We want to consider representations of the full 
Lorentz group thus we choose as the representation D 
the bispinor representation D^' ) © D^°^\ so L = 7 
in this case. Explicitly, if A E SX(2, C) and A (A) is 
an image of A in the canonical homomorphism of the 
SL(2, C) group onto the Lorentz group, we take the chiral 
form of D^'°) ® D^°'^\ namely 

D(A(A)) = { J y ^j . (15) 

The canonical homomorphism between the group 
«SX(2, C) (universal covering of the proper ortochronous 
Lorentz group L+) and the Lorentz group ~ 
SO(l, 3)o |a is defined as follows: With every four- vector 
fc M we associate a two-dimensional hermitian matrix k 
such that 

k = fc"o- M , (16) 

where <r,, i = 1,2,3, are the standard Pauli matrices 
and do = I. In the space of two-dimensional hermitian 
matrices (|16fl the Lorentz group action is given by k' = 
AkA*, where A denotes the element of the SL(2, C) group 
corresponding to the Lorentz transformation A (A) which 
converts the four-vector k to k' (i.e., fc /AI = A^ v k v ) and 

k' = 

Now, the explicit solution of the Weinberg condition 
under our choice of D (Eq. I|15|0 is given by 

+ (17, 

2 A /i + ^ V( J + ™ k mj 

V TO 

where k is given by Eq. i|16[l and k p — (k p ) fJ/ a fi with 
k p = (k°, — k). As is well known, the intertwining matrix 
v(k) fulfills the Dirac equation 

(kj - ml)v{k) = 0, (18) 

where 7^ are Dirac matrices, k-f = fc M 7 M . The explicit 
representation of Dirac matrices used in the present pa- 
per is summarized in Appendix A. 

Now we discuss the general structure of the reduced 
density matrix for s — h. We show that this matrix 
contains information about both average polarization as 
well as kinematical degrees of freedom. Recall that the 
polarization of the relativistic particle is determined by 
the Pauli-Lubanski four-vector 

W = \e^ aX P u J aXl (19) 

where P v is a four-momentum operator and J a \ de- 
notes generators of the Lorentz group, i.e., U(A) = 
expioj^ Jfj, v . We will also use the spin tensor S^ v de- 
fined by the formula || 

= -X>e tlv< , T P°W T . (20) 

711 Z 
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Now, the 4x4 reduced spin density matrix can be 
written as the following combination 

6 = i(a 7 ° + ^7 5 7° + ^7^7° + ^tVt° 

+ 2 V 4b"\ 7 "] 7 ). (21) 

Real coefficients a, 6, u^, w p , can be determined 
by calculating corresponding traces. Thus, after some 
algebra, using Eqs. ® [121131201 and (|BfHB2jl we get 



a 


= Tr(fi) = 1, 


(22) 


b 


= Tr(ifi 75 ) = 


(23) 




= Tr(n 7 ^) = i(P /1 }., 


(24) 




= |Tr(n 7M75 ) = (^)., 


(25) 




= Tr(f7|[ 7M , 7l ,]) = <V)p, 


(26) 



where (^4)^ denotes the mean value of the observable A 
in the state described by the complete density matrix p, 
(A)- = Tr (pA). Notice that the above relations are not 
accidental, since 7 ° 7 M is a canonical four-velocity opera- 
tor for the Dirac particle and |[ 7 '\ 7 I/ ] are Lorentz group 
generators in the bispinor representation. Thus, finally, 
the matrix CI = 7 ° has the following form: 



+ |(V)a^7l- (27) 



It can be also checked that in the nonrelativistic limit 
we have 



(Wo), - 0, 
(%>/, - 0, 



(28a) 
(28b) 
(28c) 
(28d) 



The formalism we have introduced above can be 
straightforward generalized to the multiparticle case. As 
an example we shall discuss briefly the reduced spin 
density matrix for two massive particles. Two-particle 
Hilbert space is spanned by vectors \a, k) ® |/3,p), where 
\a,k) is defined by Eq. (JjJJ. Therefore we define the two- 
particle unnormalized reduced density matrix as follows: 



0a'P', a p = J dfi(k) dfi(p) (a, k\ <g) (J3,p\ p \a\ k) ® \0',p) , 

(29) 

where p denotes the complete two-particle density ma- 
trix. It is obvious that the matrix l|29|) is Hermitian, 
positive-semidefinite and can be easily normalized simi- 
larly like in the one-particle case. Moreover, in the case 
of two spin 1/2 particles we define 



ci = 6>( 7 ° ® 7 °). 



(30) 



III. PARTICLE WITH A SHARP MOMENTUM 

Now let us discuss the case of the particle with a sharp 
momentum, say q, and polarization determined by the 
Bloch vector £, |£| < 1, i.e., we assume that the complete 
density matrix has the following matrix elements 

(k,m,s,T\p\p,m,s,X) 

= ^5 3 (k-q)5 3 (p-q)i(7-^- CT ) TA . (31) 

Of course the normalization factor jr^yj should be un- 
derstood as the result of the proper regularization pro- 
cedure. Now, using Eqs. (JSJ) and (|lj2fl we can find the 
corresponding matrix Ct. We have 



Cl 



4 \ m / V m / 



(32) 



where the four-vector = (W*) - p is given in this case 
by 



w 



q£ 



q(q ■ g 

q° + m 



(33) 



i.e., w is obtained from (0, m|) by applying the Lorentz 
boost L q . It should also be noted that w^q^ = 0. The 
matrix (|32() is known in the literature as the spin density 
matrix for Dirac particle 9] . 

Now, to connect the density matrix introduced above 
with some macroscopic experiments like the Stern- 
Gerlach one let us consider a charged particle with sharp 
momentum moving in the external electromagnetic field. 
We assume that the giromagnetic ratio g = 2. The mo- 
mentum and polarization of such a particle vary in time, 
thus they can be regarded as functions of its proper time 
t: 



q = q(r), € = €(r), 



(34) 



The expectation value of the operators representing the 
spin and the momentum will necessarily follow the same 
time dependence as one would obtain from the classical 
equations of motion 0, : 



dqf 



C 



■ P q + m 2 q W ° vt 13 

C -F% [F^w v + ± q e( w °F aa q a )} , (35) 

7T7 ^ L 771 Z J 



— =ZlFW + —q»{ W °F <7a q<*)].. 



(36) 



where e denotes the charge of the particle, m its mass, 
£ is the proportionality constant between the magnetic 
moment of the particle p a and w a i.e. p a = ^w°|2ll| and 
F^ u is the tensor of the external electromagnetic field, 
F a p — 'Ffj, v . Eq. lj56*|) describes Thomas precession 

of the spin vector in the electromagnetic field [1C| while 
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Eq. I|35|) allows one to determine the trajectory of the 
spinning particle moving in the electromagnetic field F^ v . 
The slow motion limit of the above equations takes the 
well-known form [TTj 



§=* q xB + W 

at m 2 



dt 



= CI x B, 



(37) 
(38) 



where we assumed that the electric component of the 
electromagnetic field is equal to zero. Eqs. (|37H38|) de- 
scribe forces acting on the particle in the Stern-Gerlach 
experiment, therefore we can really identify £ with the 
polarization of the particle. 

In this simple case of the monochromatic particle we 
can also calculate explicitly the von Nuemann entropy of 
the reduced density matrix. The matrix Q in the rest 
frame of the particle can be written as 



On — i 



1 1 
2 U 1 



(39) 



To calculate entropy we have to use the normalized den- 
sity matrix #0j but in this particular case 9q = On. Thus 
the von Neumann entropy of the state l|39|) is equal to 



Sn. 



1*1)1* 



1 



(1-1*1) In 



1-1*1 



(40) 

Now, to find the entropy in the arbitrary Lorentz frame 
we apply to the matrix 9q the Lorentz transformation 
with D(A) given by Ijl5|l and we find that entropy 
of the corresponding reduced density matrix 9' is given 
by (|4T1|l too, i.e., S§ = Sg,. Therefore for a particle with 
the sharp momentum the entropy of the reduced density 
matrix does not change under Lorentz transformations. 
However, in the case of an arbitrary momentum distribu- 
tion, the entropy of the reduced density matrix 9 is not 
in general Lorentz-invariant. 



IV. SPIN OPERATOR 

In the next section we will use our formalism to calcu- 
late the Einstein-Podolsky-Rosen-Bohm (EPR-Bohm) 
correlation function. Thus we have to introduce the spin 
operator for a relativistic massive particle. The choice is 
not obvious since in the discussion of relativistic EPR 
Bohm experiments various spin operators have been used 
[l3L fill flEl ITti IrH ITU IT9l | . However our previous consid- 
erations (Eqs. H33fl H38l) ~) as well as the classical definition 
of the relativistic spin y| suggest that the best candidate 
for the spin operator is 



S = — I W 



W°- 



(41) 



in quantum field theory |20j . It fulfills the following stan- 
dard commutation relations: 



[P»,S J ] = 0, 



(42a) 
(42b) 
(42c) 



which should be satisfied for the spin operator. Here J 1 = 
\£ijkJ^ k and one can show that it is the only operator 
which is a linear function of and fulfills relations l|42|) 

Therefore the operator corresponding to the spin pro- 
jection along arbitrary direction n (n 2 = 1) in the repre- 
sentation of gamma matrices (|A1() reads explicitly 



n . s = — I P° 

2m 



n a 
n • er 



(n x P) 



n P 

P° + 77 





(n x P) ■ er 

P er 
P er 



(43) 



where we have used Eqs. <|B6|1 . 

Eq. (|42bp implies that eigenvalues of the operator n • S 
are integers or half-integers. As one can easily check by 
direct calculation the eigenvalues of the operator l|43|l are 
equal to ±i. This observation supports our choice of the 
operator S as the spin operator. 

Now we want to express the average of the spin opera- 
tor (|41() in terms of the reduced matrix fl. One can check 
that in an arbitrary state p 



(P° + m)S) _ = - Tr (0 77 5 (J + 7 )) . (44) 



Thus a reasonable choice for the normalized average of 
the spin component is 



S = 



(P° + m)S 
(P° + m) 



Tr (Q77 5 (J + 7°)) 
2Tr + 7 )) 



(45) 



When p{k) describes a particle with a sharp momentum 
k the normalized average is simply the average of S, i.e., 
inserting reduced density matrix £1 132|) into (|45|l we get 



p(k) 



i 

2' 



(46) 



It should also be noted that in the nonrelativistic limit 
we recover the result (|46|) for an arbitrary state p 



which corresponds to the classical polarization vector £ 
(precisely to £/2) in Eq. I|33|l . This operator is also used 



(47) 
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V. QUANTUM CORRELATIONS 

Using the formalism introduced above, we now calcu- 
late the correlation between measurements of spin com- 
ponents performed by two observers, A and B, along two 
arbitrary directions, a and b, respectively. We consider 
the simplest situation in which both observers are at rest 
with respect to a certain inertial frame of reference O. 
We assume also that both measurements are performed 
simultaneously in the frame O. 

We calculate the EPR-Bohm correlation function in 
the pure state of two particles with sharp momenta 



IV) =^2c a0 \a,k) ® \p,p) . 

a/3 



(48) 



The corresponding reduced density matrix (|30fl has the 
following form: 




FIG. 1: Correlation function in the case when a _L b, k 
and a || k. b || p or a ± k, b _L p (Eq. 



lL aj3,a'/3' 



4fcW 3 (0)) 2 

v(k)v(k)j C*j°(v( P )v(p)y 



a/3 



(v(k)v(k)) T Cv(p)v(p) 



a'/3' 



(49) 



where the matrix C = (c a/ g) determines the state l|48f) 
while v(k)v(k) and v(p)v(p) are given by l|Blb|l . 

Observers A and B use observables 2a • S ® I and I ® 
2b ■ S, respectively (a 2 = b 2 = 1). Thus the correlation 
function has the form (see Eqs. 1*4*4*1) and |B8|l) 



C(a,b) 



(P° + m)a • S ® (P° + m)b • S 







- to) ® 


(p° + 








Tr 




77 5 (J 


+ 7°)) 


® (b • 


77 5 (7 


+ 7°)))" 


Tr 




f° + I) 


<8> (7 


+ o); 





{V| a • S « 



ib-S|<0) 



After some algebra we find that 



(50) 



C(a,b) = 

Tr { (b • S(p)v(p)v(p)j°) C7t (a ■ S(k)v(k)v(k)j°) T c} 



Tr 



{(v(p)v(p)j°) Ct (i;(fc)u(fc) 7 0) T C} 



(51) 



where matrices a - S(fc) and b • S(p) have the same form 
as l|43|) with n, P equal to a, k and b, p respectively. 
Now, for the sake of simplicity, we specify the state \tp). 
We choose 



C = a 



<J 2 
a 2 



(52) 



where a is a normalization constant. This choice is rather 
natural because the state described by Eqs. (p)2*|l has 
the same form for all inertial observers, namely 

17(A) ® (7(A) |V) = c «/3 1". Afc ) ® I/ 3 ' A P) . (53) 

a/3 

where we used Eq. (|15fl . Moreover in the center of mass 
frame it is an ordinary singlet state. 

Now, provided that C is given by Eq. 1 (5*21 . after 
straightforward calculation we arrive at 



C(a,b) = -a b 



(kxp) 

m 2 + kp 



(a x b) 



(a-k)(b x p) - (b-p)(a x k) 

(fc° + m)(p a + to) 



(54) 



We see that the correction to the nonrelativistic corre- 
lation function AC = C(a, b) — C n0 nrci = C(a, b) + ab is 
of order (3 2 , where j3 = -|, v denotes the velocity of the 
particle, and c the velocity of light. Let us note first that 
when momenta of both particles are parallel or antipar- 
allcl the correlation function has the same form as in the 
nonrelativistic case. This result differs from Czachor's 
results [2^ H3- The reason is that we use a different, 
and in our opinion more adequate, spin operator. 

Now let us consider the configuration in which the non- 
relativistic correlation vanishes (alb). For simplicity 
let us assume also that |k| = |p| and a || k, b || p or 
a 1 k, b 1 p. In such fixed configurations the correla- 
tion function has the very simple form 



C(a,b) = AC = 4-4 

Pq + TO 



I? 



2-I3 2 



(55) 



Dependence of the above correlation on (3 is depicted in 
Fig. HQ Notice that l|55|) was also obtained by Czachor 
|14| but for a different configuration. 
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VI. CONCLUSIONS 

To conclude, we have constructed a Lorentz-covariant 
reduced spin density matrix for a single massive particle. 
It contains not only information about average polariza- 
tion of the particle but also information about its average 
kinematical state. We have also showed that this matrix 
has the proper nonrelativistic limit. 

Our results shows that we can define a Lorentz- 
covariant finite-dimensional matrix describing polariza- 
tion of a massive particle. However in the relativis- 
tic case (contrary to the nonrelativistic one) we cannot 
completely separate kinematical degrees of freedom if we 
want to construct a finite-dimensional covariant descrip- 
tion of the polarization degrees of freedom. 

With help of our covariant formalism we have also cal- 
culated the correlation function in the EPR-Bohm type 
experiment with massive relativistic particles. We have 
showed that relativistic correction AC to the nonrelativis- 
tic correlation function C nonre ] = —a • b vanishes when 
momenta of both particles are parallel or antiparallel, 
i.e., in the standard configuration of EPR-Bohm type 
experiments. We have found also the configurations in 
which the nonrelativistic correlation vanishes while the 
relativistic correction AC survives and is of order /3 2 (Eq. 
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APPENDIX A: DIRAC MATRICES 

In this paper we use the following conventions. Dirac 
matrices fulfills the condition 7^7 l/ + -y^y = 2g^ v where 
g^ v = diag(l, —1, — 1, —1) denotes Minkowski metric ten- 
sor; moreover we adopt the convention e 0123 = 1. We use 
the following explicit representation of gamma matrices: 

^-(?0- '-('7). 

(Al) 



where a = (<ti, 02, 03) and Oi are standard Pauli matri- 
ces. 

APPENDIX B: USEFUL FORMULAS 

The matrix (|17|l is normalized as follows 

v(k)v(k) = I, (Bla) 
v(k)v(k) = ^(k7+ml), (Bib) 

where v(k) = v^(k)^°. Moreover it can be verified that 
it fulfills the following relation 

vtkWvik) = —I. (B2) 
m 

Vectors \a, k) fulfill the orthogonality relation: 

(a, k\p,p) = 2fc°<5 3 (k - p) {v(k)vHk)) 0a , (B3) 
and one can check that 

I = J2jMk)^ a \a,k){P,k\, (B4) 

a/3 

(v(k)v(k)) \p,k) = \a,k). (B5) 

V / a/3 

In the representation of gamma matrices (jAlJl we have 

*-WH ft <B6b) 

It can be also checked, that when 

(P) =Tr(n/), (B7a) 
(d) =Tr(Op), (B7b) 

we have 

(f®&) =TV(fi(/® g )), (B8) 

where in Eqs. <|BT|) and (|B8|) p and f2 are complete and 
reduced density matrices for one and two particles, re- 
spectively. 



[1] A. Peres, P. F. Scudo, and D. R. Terno, Phys. Rev. Lett. 

88, 230402 (2002). 
[2] P. Caban, K. A. Smoliriski, and Z. Walczak, Phys. Rev. 

A 68, 044101 (2003). 
[3] M. Czachor, Phys. Rev. Lett. 94, 078901 (2005), quant- 

ph/0312040. 



[4] A. Peres, P. F. Scudo, and D. R. Terno, Phys. Rev. Lett. 
94, 078902 (2005). 

[5] S. Weinberg, Phys. Rev. 133, B1318 (1964). 

[6] S. Weinberg, Phys. Rev. 134, B882 (1964). 

[7] A. O. Barut and R. Raczka, Theory of Group Represen- 
tations and Applications (PWN, Warszawa, 1977). 



7 



[8] J. L. Anderson, Principles of Relativity Physics (Aca- 
demic Press, New York, 1967). 

[9] W. B. Berestetzki, E. M. Lifschitz, and L. P. Pitajew- 
ski, Relativistic Quantum Theory, vol. 1 (Nauka, Moscow, 
1968). 

[10] V. Bargmann, L. Michel, and V. L. Telegdi, Phys. Rev. 

Lett. 2, 435 (1959). 
[11] J. P. Costella and B. H. J. McKellar, Int. J. Mod. Phys. 

A 9, 461 (1994). 
[12] H. C. Corben, Phys. Rev, 121, 1833 (1961). 
[13] D. Ahn, H. J. Lee, Y. H. Moon, and S. W. Hwang, Phys. 

Rev. A 67, 012103 (2003). 
[14] M. Czachor, Phys. Rev. A 55, 72 (1997). 
[15] J. Rembielihski and K. A. Smoliriski, Phys. Rev. A 66, 

052114 (2002), quant-ph/0204155. 



[16] D. Lee and E. Chang- Young, New J. Phys. 6, 67 (2004). 
[17] H. Li and J. Du, Phys. Rev. A 68, 022108 (2003). 
[18] H. Terashima and M. Ueda, Int. J. Quant. Inf. 1, 93 
(2003). 

[19] H. Terashima and M. Ueda, Q. Inf. Comput. 3, 224 
(2003). 

[20] N. N. Bogolubov, A. A. Logunov, and I. T. Todorov, 
Introduction to Axiomatic Quantum Field Theory (W. 
A. Benjamin, Reading, Mass., 1975). 

[21] For the particle with charge e and giromagnetic ratio g 

[22] The Czachor's result can be obtained in our framework 
by calculating the appriopriately normalized average of 
a • W <g> b • W. 



